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Abstract 

We calculate one- loop renormalization factors of generic As = 2 four-quark operators for domain- wall QCD 
with the plaquette gauge action and the Iwasaki gauge action. The renormalization factors are presented in 
the modified minimal subtraction (MS) scheme with the naive dimensional regularization. As an important 
application we show how to construct the renormalization factors for the operators contributing to /■sT" — K'^ 
mixing in the supersymmctric models with the use of our results. 
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I. INTRODUCTION 
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Experimental studies of — mixing provide us an opportunity to deduce the indirect 
CP violation parameter eK- In the Standard Model (SM), the low-energy effective Hamiltonian 
contains the dimension-six four-quark operator Oll = s7/i(l — 75)^^' s^^J.{^ ~ lb)d-, and its hadronic 
matrix element {K^\Oll\K^) is required to determine ex from the experimental results of — K^ 
mixing. On the other hand, the physics beyond the SM involves the four-quark operators with 
more general chiral structures. For example, the relevant operators in supersymmetric models 
are0| 

(1) 
(2) 
(3) 
(4) 
(5) 

where a and h are color indices. 

Lattice QCD should be an ideal tool to determine the above matrix elements from the first 
principles. In the past decades, much effort have been devoted to the calculation of {K^\Oll\K^) by 
employing various quark and gauge actions 0, 0, Q, 0, 0, Q, 0, 0, Hfll . . On the contrary, the matrix 
elements for the operators of Eqs.(P) — © are less studied compared to the case of Qj- 
Our aim is at a detailed study of them with the use of the domain- wall quark formulation in lattice 
QCD, which is expected to realize full chiral symmetry at finite lattice spacing up to exponential 
fall-off of the explicit chiral symmetry breakingcontributions^Jj and have been successfully 
applied to the calculation of (i?°|OLL|K°) 0,0,13. 

The matrix elements calculated on the lattice should be converted to those defined in some 
continuum regularization scheme [e.g., the modified minimal subtraction scheme (MS)]. This is 
achieved by the finite renormalization relating the lattice composite operators to the continuum 
counterparts defined in some reguarization scheme. In this paper we present the perturbative results 
of the renormalization factors for the complete set of As = 2 four-quark operators consisting of 
physical quark fields in the domain- wall QCD (DWQCD). With the use of our results one can 
obtain the renormalization factor for an arbitrary As = 2 four-quark operator. This work is an 
extention of Eefs.jl^. Il7l|. where the renormalization factors for the four-quark operators relevant 
in the SM are evaluated. 

This paper is organized as follows. In Sec. Owe introduce the quark and gauge actions and the 
corresponding Feynman rules. Section ITTll is devoted to explain our calculational procedure of the 
renormalization factors for the complete set of As = 2 four-quark operators. In Sec. IIVI with the 
use of our results, we construct the renormalization factors for the SUSY operators in Eqs.(P) — Q. 
We briefly discuss the mean field improvement in Sec. Our conclusions are summarized in 
Sec. EH 

The physical quantities are expressed in lattice units and the lattice spacing a is suppressed 
unless necessary. We take SU(A'') gauge group with the gauge coupling g and the second Casimir 
Cp = {N"^ — 1)/(2A^), while = 3 is specified in the numerical calculations. 
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II. ACTION AND FEYNMAN RULES 



We take Shamir's domain- wall fermion action [l^ given by 



Ns 



n s=l 



2 



+mY[il;{n)NsPR'>P{n)i + V("-)i^LV'(?^)AfJ , (6) 

n 

where n is a four-dimensional space-time coordinate and s is a fifth-dimensional or "flavor" index 
bounded as 1 ^ s ^ A'^,. In this paper we conventinally take Ng ^ 00 limit to avoid complications 
arising from the finite Ng such as mixing of operators with different chiralities. The domain-wall 
height M is a parameter of the theory which we set < M < 2 to realize the massless fermion 
at tree level. Pji/l is a projection operator Pr/l = (1 i 75)/2 and the Wilson parameter is set to 
r = — 1. The "physical" quark field are defined by the boundary fermions in the fifth dimensional 
space: 

qin) = PRiP{n)i+PLtl;{n)N^, (7) 
q{n) = ^{n)N^Pn+^{n)iPL, (8) 

whose mass is given by m. Our renormalization procedure is based on the Green functions con- 
sisting of only the "physical" quark fields. 

For the gauge part of the action we employ the following form in four dimensions: 

1 



'S'gluon — 9 

5 L 



CO Yl T^^pi + ci Yl Trf/rtg + C2 Trf/ehr + C3 Y Trf/pigJ, (9) 

plaqucttc rectangle chair parallelogram 



where the first term represents the standard plaquette action and the remaining terms are six-link 
loops formed by a 1 x 2 rectangle, a bent 1x2 rectangle (chair) and a 3-dimensional parallelogram. 
The coefficients cq , . . . , C3 satisfy the normalization condition 

Co + 8ci + 16c2 + 8c3 = 1. (10) 

The RG improved gauge action is defined by setting the parameters to the value suggested by an 
approximate renormalization group analysis. In the following we adopt ci = — 0.331, C2 = C3 = 
flwasaki) (l^ for the RG improved gauge action in addition to ci = C2 = C3 = (plaquette). With 
these choices of parameters the RG improved gauge action is expected to realize smooth gauge 
field fluctuations approximating those in the continuum limit better than with the unimproved 
plaquette action. 

Weak coupling perturbation theory is developed by expressing the link variable in terms of the 
gauge potential 

[7^,^ = exp{igA^{x + ^fi)) (11) 

and expanding in terms of the gauge coupling. We adopt a covariant gauge fixing with a gauge 
parameter a defined by 

^GF = E^(^/^^m(^+^/^))'' (12) 
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where V^/„ = {fn+fi - fn)- 

The free part of the gluon action takes the form in momentum space 



-k), 



where 



Gfj,u{k) — k^ku + ^ '^ jkpS^u f^fi^pu)Q^ipkp 



with 

and Qpiy is defined as 

Qfiu = (1 - Sfiu) 



2sm — 
2 



1 - (Ci - C2 - C3){kl + kl) - (C2 + C3)A:^ 



The gluon propagator can be written as 



D^,{k) = {k^y 



kpku + ^ '^ jkqSp^u ki,6fj_(^)ko-A, 

(T 

(1 ~ ■Ap.u)kpki, + dp^y ^ ^ k^Apfj 



(l-a) 



kn,k^ 



(l-a) 



(A;2)2 
kpky 



(13) 
(14) 

(15) 
(16) 

(17) 
(18) 



where Ap^ is a function of (7^,^ and kp , whose form we refer to the original literatures 3) ^ ■ Iii this 
paper we will adopt the Feynman gauge(a = 1) without loss of generality, since the renormalization 
factors for the composite operators do not depend on the choice of the gauge fixing condition. 

Quark-gluon vertices are also identical to those in the A^^ flavor Wilson fermion. We need only 
one gluon vertex for our present calculation: 



V^^{k,p)st = -igT^{-fpCos{-kp/2 + pp/2) - irsm{-kp/2 + pp/2)}dst, 



(19) 

where k and p represent incoming momentum into the vertex (see Fig.l of Ref.[20|) and T^{A = 
1, . . . , A'^^ — 1) is a generator of color SU(A^). 

The fermion propagator originally takes Ng x Ng matrix form in s-flavor space. In the present 
one-loop calculation, however, we do not need the whole matrix elements because we consider Green 
functions consisting of the physical quark fields. The relevant fermion propagators are restricted 
to the following three types: 



{qi-p)qip)) 
{q{-p)ip{p, s)) 

{ip{-p,s)q{p)) 



— sinp^ + (1 — T4^e ")m 
-(1 - e°'W) + m2(l - We-°') " 

[i-fp sinpp - m(l - VFe"")] (e 



Sqip), 



(20) 



PR + e 



-q(s-1) 



Pl) 



m 



[i-fpsmpp-m{l - We~")]-F e~"(e-"("-^)PR + e 



~a(Ns~s) 



F 



(e 



Pl + e-^^'^-i^Pfi) [i-fp sinpp - m(l - T^e"")] 



+ (e-"{^-i)p^ + e-"(^^-^)PR)e-° \m[i-fpsmpp - m{l - We 



Pl)121) 



(22) 
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with 



W = 1 -M-r -cos p^), (23) 

cosh(a) = ^ ^, (24) 

F = l-e"T^-m2(l-M^e-"), (25) 

where the argument p in the factors a and W is suppressed. 

In the perturbative calculation of Green functions we assume that the external quark momenta 
and masses are much smaller than the lattice cutoff, so that the external quark propagators can be 
expanded in terms of them. We have the following expressions as leading term of the expansion: 

i]p + [I — WQ)m 

(#s)cxt(p) = {qq)ip){w'o~'PL+w^'-'PR), (27) 
(^sg)cxt(p) = {w^o~'PR + <'~'PL){QQ)ip), (28) 

where wq = 1 — M. 



III. RENORMALIZATION FACTORS FOR GENERIC FOUR-QUARK OPERATORS 

We consider the complete set of parity-conserving four-quark operators 0, 0]: 

^ O^y + O^^, (29) 

^ 0$v - O^A^ (30) 
0%-0$p, (31) 

0% + 0$p, (32) 

(33) 

where 

1 r. 



Qt 


— ^VV+AA 


Qt 


— ^VV-AA 


Qt 


— ^SS-PP 


Qt 


— ^SS+PP 


Qt 


= 



Otr - - 



{ifiTq^) iqFql) ± (g^Fg^) (g^Fg^) (34) 



with r = {1, 7;^, cr^jy, 7/i75, 75} = {S, V, T, A, P}. a, b denote color indices and summation over them 
is assumed. We should note that qi (1=1,2,3,4) are the physical quark fields defined by Eqs.Q and 
(jHJ. The parity-conserving parts of Eqs.(^ — @ are expressed as linear combinations of the above 
operators. We find that Ref.Q] employs another choice of basis called Fierz eigenbasis, which is 
more convenient in considering the Fierz transformation of the four-quark operators. 
We consider the following Green functions: 

{Qt)%% ^ {Qf{<lim2%{q,f,{qA)\), (35) 

where a, /3, 7,(5 are spinor indices and i,j,k,l are color indices. It is convenient to decompose the 
above Green functions as 

{OtXH.s ^ (Orr('?i)i('?2)^^(g3)^(94)^) 

= \mrq^2){'i^ql) ± (^-?rg2)(^1rg^)](9i)L(92)^^(g3)^(^4)^> 

= I [((-zOL^-^rgf (g-2)^^(g3)^g1rg^(Q-4)i) T (('7i)L9?rg2(g-4)kg3)^9lrg^(g2)p] . (36) 
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After trancating the external quark propagators from {Qf) with the multiplication oii]f+{l—WQ)m, 
we obtain the vertex functions, which is written in the following form up to the one-loop level: 



2n4/a{0)± . .{l)±xii;fcZ 



(37) 



where the superscript (n) refers to the nth loop level and the subscript i identifies the operators 
(|29|) — (IHI-^j) ■ The trivial factor (1 — Wq)^, which originates from the external quark propagators, is 
factored out for convenience. Since the renormalization factor does not depend on the external 
momenta pi, we suppress them. 



The tree level vertex functions A-''^^ are given by 

i = 1, ^[V(^V + A^A^{VQV + AQA)]o,p,^s[l^ir''''\ 
- 2, l^[V 0V - A0A^{V QV - AQA)]a,/3--ys[l^lf'''^, 



[S0S-P^Pt{SQS-PQ P)Uf3,^s[l^l] 

[S0S + P0Pt{SQS + PQ P)]a/3;75 [1(8)1] 



ij;kl 



ij:kl 



[T®TTTQT]a(i;^s[l^l] 



ij;kl 



(38) 
(39) 
(40) 
(41) 
(42) 



where (8), act on the Dirac spinor space as [F T]a/3-'yS = {T)a/3{T).yS, [T T]af3-'y5 = (r)c5(r)^/3, 
and 0,0 on the color space as [101]*-'-'""' = 6ij6ki, [101]*-'''^' = 6ii5kj- 

Now let us consider the one-loop vertex corrections depicted by six diagrams in FigQ Their 
total contribution yields the vertex function at one-loop level: 



(i)± 



A 



A 



(i)± 



A 



(i)± 



A 



(i)± 



^ d^k 
^ d^k 
^ d'^k 



{lyv + Iaa) + ^ . . . , + [lyv + ^Aa) 
{^VV ~ -^1a) + ) • • • ^+{^VV ~ ^aa) 

i^ss ~ Ipp)+^ . . . , + {iss ~ Ipp) 

(/^5 + /^p)+,...,+(/^'s + /|,'p) 

Itt~^i . . . , +Itt 



(43) 
(44) 
(45) 
(46) 
(47) 
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where 

Irr = \g^{T^T^®l)K [{^alpT^pla ® r(^yl„^ sin^ sin^ kp + sin^ ^ cos^ ^ sin^ kp) 

(7 

+ r0r(4Ai-^J]Af,i + ^J]sin^fc„+T)}T{®^0}], (48) 
-^fr = (r^®r^)K ^{7Q7/3r <8> r7^7„ (^-4a/3 sin^ A;„ sin^ % + X] "^"/^ ^^"^^ Y '^^^^ Y '^^^ ^Z^) 

cr 

+ r0r(4A2-l^Af,i + ij]sin4A;.+T)}T{®-0}], (49) 

Of cr 

-^rr = -\9^{T^®T^)K [{7a7/3r 7Q7/3r(^A«/3 sin^ k^, sin^ /c/3 + ^ sin^ cos^ -y sin^ kp) 

a 

+ r0r(4A2-l^A?,i + iXsin2fc. + r)}T{0--0}] (50) 

(T cr 

with 





4 1 


(51) 




= - Wo, 


(52) 


F 


= e"" - W, 


(53) 


T 




(54) 


As 




(55) 


Ai 




(56) 


AM 

^1,1 




(57) 



Summation over repeated indices is assumed. The above expressions are obtained with the aid of 
useful formula presented in Refs. , [itI . Eo| . We should note that the other three contributions 
ja^b ,c Figs, la', lb', Id are equal to 1^^'^ respectively. After all, the total contribution 
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becomes 



{^ss ~ Ipp)+-> • • • 1 

(/gs + /^p) +,..., +(/S'5 + /|>'p) 



xK[{^ T l)AsP + {2Cf - ^ ± l)AvA + 2CfT], (58) 
.+{lvv -Iaa) = 9^[V®V-A®At{VqV-Aq A)]^p.,-^s[l®lV-'"' 
xK[-^AsP + {2Cf + ^)AvA + 2CfT] 
+g'^[S®S-P®PTiSQS-PQ P)Uf3:^s[l®iy^--''^ 
xK[t2Asp±2Ava], (59) 

+{iss - ipp) = As 0S-P0Pt{sqs-pq p)U.,^s[mr--''' 

xK[2CfAsp + 2CfT], (60) 
g^[S(^S + P^PT{SQS + PQ P)U:^s[l^ir'-''^ 
xK[{2Cf T l)AsP ± AvA + 2CfT] 

(61) 



xi^[(^T^)^.P + (-^±^)AvA)], 



ij:kl 



xK[i-^ T l)AsP + (^ ± l)AvA + 2CfT] 
+g'^[S®S + P(g)PTiSQS + PQ P)]o,p.,^s[l®l] 
xK{^±1)[Asp-Ava], 



ij:kl 



where we use the Fierz rearrangements: 

1 



[S S] [101] 
[V V] [101] 

[T r] [101] 

[^ A] [101] 
[P P] [101] 



[5 05 + y0y-T0r-^0^ + P0P] [101] 

[45 5 - 2y y - ^ - 4P p] [101] , 

[-65 5 - 2T r - 6P P] [101] , 

[-45 5 - 2V y - A + 4P P] [101] , 

[5 05-y0y-T0r + ^0^ + P0P] [101] 



(62) 

(63) 
(64) 
(65) 
(66) 
(67) 



Note that these formula do not include Fermi statistics. 
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Comparing the one-loop results to the tree level ones we obtain 



I + 9^^^^{t{Asp) + {N± 2){Ava) + {N± 1)(T)} 



± 



A 



A3^ 



Ar 



N 

1 + g^{-j^{Asp) + {2Cf + j^){AvA) + 2Cf{T)} 
+g^{T2{Asp)±2{AvA)}Ar^, 



A 



{0)± 

1 ' 



A 



(0)± 



(0)± 



l + g^2CF{{Asp) + {T)} A 



1 + g^^m^ tN- l){Asp) ± N{Ava) + (N^ - 1)(T)} 



+9 



N 
3N 



(0)± 



{{Asp) - {Ava)}aP^, 



1 + g'ii-^ T l){Asp) + ± + 2Cf{T)} 



A 



+<7^(-±l){(A5p)-(^yA)}A, 



(0)± 



with 



K{k)X{k) 



(68) 

(69) 
(70) 

(71) 

(72) 

(73) 



for X = T, AyA, Asp- We remark that Cf{T + ^y^) and Cf{T + Asp) correspond to the one- 
loop vertex corrections to the (axial) vector and the (pseudo) scalar density. Ay a and Asp are 
expressed as follows: 



AvA = 4A2-A^i + 4.Xo> 
Asp = lOAgA^n 



with 



■ 2 

sm — . 
2 2 



(74) 
(75) 



(76) 



The expressions of Eas. H68p and ()7U() show important properties of Qi and Q3 in DWQCD for- 
malism: the one-loop vertex corrections are multiplicative. This is contrary to the Wilson case, in 
which the mixing operators with different chiralities appears at the one-loop level. The results for 
Qf in Ea. H68() are already obtained in Ref. jl^ . 

Combining the contribution from the quark self-energy evaluated in Ref. [i^l and the vertex 
corrections, we obtain the lattice renormalization factors: 



where 



Z2 



1 + 
1 - 
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Ci.[log(Aa)2 + Si], 
S3, 



167r2 
2wo g^CF 



\ — w\ IGvr^ 



± 



9 



167r2 L 



[75log(Aa)2 + 



(77) 

(78) 
(79) 

(80) 
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A is the fictitious gluon mass introduced to regularize the infrared divergences. Si and S3 are 
finite parts of the renormalization factors of quark wave function and overall factor (I — Wq), whose 
numerical value is given in Ref. 17 1. The matrix is expressed as 

/vf^ \ 



'11 

^22 ^23 

vf^ 

vt 





















^t5 



V 



(81) 



4 vfj 



whose components are 
+ _ 167r2(Ar-l) 



'11 



'11 



vt2 



^23 



^33 



^44 



^45 



■"55 



^ [-{{Asp)) + {N + 2){{Ava)) + {N + 1)(T)] +7+ logvr^ 
167r2(iV + 1) ^^^^^^^^ ^ _ 2)((^^^)) + _ 1)^^)] log ^2^ 



N 



[-{{Asp)) + N\{Ava)) + (N^ - l){T)]+ji,log7T' 



levr^ 

Tl67r2[2((A5p)) - 2((^y^))]+7|log7r2 
167r2(iV + l)(iV - 1) 



[{{Asp)) + {T)]+^i,log7r' 



lOvr^ 



[((Asp))-((^yA))]+741logv^^ 



167r2(2 T A^) 



3N 



(82) 

(83) 

(84) 
(85) 

(86) 
(87) 
(88) 



-^[ ^ ((^5P)) + ^ ((^va)) + (iV' - l)(r)]+75± log7r2,(89) 

167r2(2 ±iV) 



N 



[{{Asp)) - {{AvA))]+7tlog7r' 



(90) 



with 



/7n 



\ 




72I 72I 

7^3 

7± 7± 

V 75I 751/ 



/ 


{NTl)iNT2) 
N 











\ 







{N+2)[N~2) 
N 


T6 
















(7V+1)(7V-1) 

^ AT 



















4N'^^3N-4: 
N 


2^N 
N 


V 











6±3N 
N 


/ 



(91) 



The infrared singularity in {Ax) is subtracted as 



{{Ax)) 

where csp = 4: and cyA = 1- 



. (2vr)^ 



i^(fc)^x(fc)-cx 



0(^2 -A;2) 
(A;2)2 



(92) 
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The lattice operators and the continuum ones defined in the MS scheme with the naive dimen- 
sional reguralization (NDR) are related as 



with 



(1 



(93) 



EfS(NDR) 



(NDR) 



±MS 



^ 167r2 L 



(7ij - 2CFSij) log{naf + 



„±MS 



MS 



,+MS 

hi 



(NDR) = -5, 



V. 



±MS „,±MS 



22 



23 



±MS „,±MS 



V 



32 






33 












V 





_ q_ 

±MS „,±MS 



44 



45 



,±MS „,±MS 



-"54 



55 



7 



-MS 



(NDR) = 6 



=^2 3 








V 






3 

1 



17 






=Fl 
^2/ 



(94) 
(95) 
(96) 

(97) 



where ^g^^ are for the case of = 3. The continuum counterparts of the wave- function renor- 
malization factor and the vertex corrections are evaluated by employing the same gauge fixing 
condition and the same infrared regulator as the lattice case. Here it should be noted that the 
one-loop vertex corrections require to specify the evanescent operators, which originates from the 
property that the Fierz transformation cannot be defined in the NDR scheme jl^]. We employ the 
following evanescent operators: 



E- 



NDR 
VV 

^SS 
Epp 

jrnNDR 



a7/3 



-{2(1- 

-{(4- 
{6(2 



-2D + 2)V ®V + 2(1 - D)A ® A}, 
D)V + {D^ - 2D + 2)A ® A}, 
- 5L> + 8)5 S* + (4 - D)P (g) P + (2 



D)S^S + {D^ - 
- D)S ^ S + 6{2 



5D + 
■D)Pi 



)P(E)P + {2 



D)Ti 
- D)Ti 



(98) 

(99) 

.T}, (100) 
DT}, (101) 



P + {D'^ -2D + A)T (g) r},(102) 



where D is the reduced space-time dimension. 

Numerical values for Vij and Zij are evaluated by momentum integration, which is performed 
by a mode sum for a periodic box of a size after transforming the momentum variable through 
kfj, = — sinq'^. We employ the size L = 64 for integrals, and numerical error is estimated by 
varying L from 64 to 60. The results are presented in Table I,. . . ,VI for the plaquette gauge 
action and in Table VII,. . . ,XII for the Iwasaki gauge action as a function of M. These values are 
evaluated with = 3. 
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IV. RENORMALIZATION FACTORS FOR SUSY OPERATORS 



Let us consider the parity-conserving parts of the SUSY operators in Eqs.(^ — which are 
relevant for — mixing. They are related to the operators of Eas.(| ^ — as 

01 = Qt, (103) 

02 = Qt, (104) 

03 = -^iQt-Qt)^ (105) 

04 = Qt, (106) 

05 = -^Qt, (107) 

where we take qi = q-i = s and q2 = q^ = d in Eq. plj) . With the use of these relations, the 
renormalization factors for Oi are expressed as 

= n — \w^[zti{HO'r{i/a)], (108) 

^^(/^) = n W2 [(^4'4(^«) + Zt,{lia))0't{\la) + 2Z+ (HOf (1/a)], (109) 

^Til^) = J. \^\{Zt^{iia)-Zt^{i,a))Of(\la) 

-\{Zt,{iia) + Z+ (/.a) - Z5+ (/.a) - Z5+ (^a))Oi-'(l/a)], (110) 

^F^(^) = n \^\Ztz{P^'^)OT(M^) - 2Zt2{m)Of(\la)\ (111) 

= n ^T2^[^2'2(HOf (l/«)-^^2'3(M«)C^f (l/«)], (112) 

where Zf- are given in Eg . (1941) . 



V. MEAN FIELD IMPROVEMENT 

Let us briefly explain the mean-field improvement on the renormalization factors of the four- 



quark operators in Eas. H29() — (|33() . We follow the discussion in Sec. VI of Ref.|l7t|. The renormal 



ization factors with the mean field improvement for the four-quark operators are given by 



= 7r^4272^S(/^«)^f''(V«), (113) 

(l-Wgj^Z^ 



zt = + ' 



2 



167r2 



(7i,-2Ci.5i,)log(H'+4''^ , (114) 



^±MF ^ ^±MS _ ^± + 5,,2CF{^f^ - Si) + S^.Wtt'CpTmf, (115) 
where Tmf is the one-loop correction to the mean field factor defined by 

n = l-g'CF^ + ... (116) 

and Wo = I — M with M = M — 4(1 — u). Note that the overall factor (1 — uJq)'^Z^ is mean-field 
improved following the description in Ref . [l7| . 
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Now it is instructive to evaluate the magnitude of the renormahzation factors in current repre- 
sentative simulations. We take Ref.|23] as an example, where a quenched simulation of domain- wall 
QCD was made at /5 = 6.0 for the plaquette gauge action with M = 1.8 and /? = 2.6 for the Iwasaki 
gauge action with M = 1.8. The mean field improved MS coupling 9|jg(l/o^) the scale /x = 1/a 
is obtained from 

-^(l/a)=P^ + dg + Cp (117) 



for the plaquette gauge action, while 

,/3 



— (1/a) = (coP + 8cii?)^ + + CO • Cp + 8ci • cri (118) 



for the Iwasaki gauge action, where the values of dg and Cp are listed in Table XVI of Ref.|17l|. P 
denotes the expectation value of the plaquette and R for the 1x2 rectangular. Their values are 
taken from Ref.0|. The domain-wall height M = 1.8 is replaced with M = 1.3112 (plaquette) and 
M = 1.4198 (Iwasaki), respectively, according to M = M — 4(1 — n) with u = P^/^. The mean field 
improved values for the renormahzation factors Zf^ at the scale fi = 1/a are given in Table IXlTll 
We find reasonable magnitude of corrections to the tree-level results for all the renormahzation 
factors. 



VI. CONCLUSION 



In this paper we have evaluated the renormahzation factors for the generic four-quark operators 
at the one-loop level in DWQCD with the plaquette and Iwasaki gauge actions. As an application 
we explain how to construct the renormahzation factors for the SUSY operators by using our 
results. We also show that taking the parameters employed in Ref.ji^l the numerical values for the 
renormahzation factors result in reasonable magnitude with the mean-field improvement. 
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TABLE I: Numerical values for Vij (i, j = 1, . . . , 5) as a function of M with plaquette gauge action. 



M 


<i 




^'11 






"23 


^^23 


0.1 


13.5696(73) 


11 


5369(73) 


13 


2308(18) 


2.0327(66) 


-2.0327(66) 


0.2 


13.0548(73) 


12 


5870(73) 


12 


9768(18) 


0.4678(66) 


-0.4678(66) 


0.3 


12.6404(73) 


13 


4380(73) 


12 


7733(18) 


-0.7976(66) 


0.7976(66) 


0.4 


12.2775(72) 


14 


1882(72) 


12 


5960(18) 


-1.9107(65) 


1.9107(65) 


0.5 


11.9450(73) 


14 


8801(73) 


12 


4342(18) 


-2.9352(66) 


2.9352(66) 


0.6 


11.6307(72) 


15 


5380(72) 


12 


2819(18) 


-3.9073(65) 


3.9073(65) 


0.7 


11.3269(73) 


16 


1781(73) 


12 


1354(18) 


-4.8512(66) 


4.8512(66) 


0.8 


11.0275(72) 


16 


8129(72) 


11 


9918(18) 


-5.7853(65) 


5.7853(65) 


0.9 


10.7276(74) 


17 


4523(74) 


11 


8484(18) 


-6.7247(66) 


6.7247(66) 


1.0 


10.4225(73) 


18 


1071(73) 


11 


7033(18) 


-7.6845(65) 


7.6845(65) 


1.1 


10.1073(73) 


18 


7869(73) 


11 


5539(18) 


-8.6796(66) 


8.6796(66) 


1.2 


9.7767(75) 


19 


5039(75) 


11 


3979(19) 


-9.7271(67) 


9.7271(67) 


1.3 


9.4243(73) 


20 


2719(73) 


11 


2322(18) 


-10.8475(66) 


10.8475(66) 


1.4 


9.0419(75) 


21 


1082(75) 


11 


0529(19) 


-12.0664(67) 


12.0664(67) 


1.5 


8.6183(73) 


22 


0373(73) 


10 


8548(18) 


-13.4191(65) 


13.4191(65) 


1.6 


8.1375(74) 


23 


0928(74) 


10 


6301(18) 


-14.9553(66) 


14.9553(66) 


1.7 


7.5747(73) 


24 


3277(73) 


10 


3669(18) 


-16.7531(66) 


16.7531(66) 


1.8 


6.8864(73) 


25 


8324(73) 


10 


0440(18) 


-18.9460(66) 


18.9460(66) 


1.9 


5.9812(73) 


27 


7944(73) 


9.6168(18) 


-21.8132(66) 


21.8132(66) 



TABLE IL Numerical values for Vij (i, j = 1, . . . , 5) as a function of M with plaquette gauge action. 



M 




'^3*3 






«44 


^'45 


"45 


0.1 


10 


182(12) 


11 


1981(84) 


9.165(15) 


0.11293(37) 


-0.5646(18) 


0.2 


12 


275(12) 


12 


5090(84) 


12.041(15) 


0.02599(37) 


-0.1299(18) 


0.3 


13 


970(12) 


13 


5709(84) 


14.369(15) 


-0.04431(37) 


0.2216(18) 


0.4 


15 


462(12) 


14 


5067(83) 


16.417(15) 


-0.10615(36) 


0.5307(18) 


0.5 


16 


837(12) 


15 


3693(84) 


18.305(15) 


-0.16307(36) 


0.8153(18) 


0.6 


18 


143(12) 


16 


1893(83) 


20.097(15) 


-0.21707(36) 


1.0854(18) 


0.7 


19 


412(12) 


16 


9866(84) 


21.838(15) 


-0.26951(37) 


1.3476(18) 


0.8 


20 


670(11) 


17 


7771(83) 


23.562(15) 


-0.32141(36) 


1.6070(18) 


0.9 


21 


935(12) 


18 


5731(85) 


25.298(15) 


-0.37359(37) 


1.8680(18) 


1.0 


23 


230(12) 


19 


3878(84) 


27.072(15) 


-0.42692(36) 


2.1346(18) 


1.1 


24 


573(12) 


20 


2335(84) 


28.913(15) 


-0.48220(37) 


2.4110(18) 


1.2 


25 


989(12) 


21 


1250(86) 


30.852(15) 


-0.54040(37) 


2.7020(19) 


1.3 


27 


504(12) 


22 


0798(84) 


32.927(15) 


-0.60264(37) 


3.0132(18) 


1.4 


29 


153(12) 


23 


1193(86) 


35.186(15) 


-0.67036(37) 


3.3518(19) 


1.5 


30 


983(12) 


24 


2739(84) 


37.693(15) 


-0.74550(36) 


3.7275(18) 


1.6 


33 


063(12) 


25 


5853(85) 


40.541(15) 


-0.83085(37) 


4.1542(18) 


1.7 


35 


496(12) 


27 


1199(84) 


43.873(15) 


-0.93073(37) 


4.6536(18) 


1.8 


38 


463(12) 


28 


9900(84) 


47.936(15) 


-1.05255(37) 


5.2628(18) 


1.9 


42 


337(12) 


31 


4299(84) 


53.243(15) 


-1.21184(37) 


6.0592(18) 
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TABLE III: Numerical values for (i, j = 1, . . . , 5) as a function of M with plaquettc gauge action. 



M 






"55 


^^54 


«54 


0.1 


14.8119(33) 


12 


7791(33) 


-1.6939(55) 


0.3388(11) 


0.2 


13.3406(33) 


12 


8729(33) 


-0.3898(55) 


0.0780(11) 


0.3 


12.1530(33) 


12 


9506(33) 


0.6647(55) 


-0.1329(11) 


0.4 


11.1099(32) 


13 


0206(33) 


1.5922(54) 


-0.3184(11) 


0.5 


10.1512(33) 


13 


0864(33) 


2.4460(55) 


-0.4892(11) 


0.6 


9.2429(33) 


13 


1502(33) 


3.2561(54) 


-0.6512(11) 


0.7 


8.3622(33) 


13 


2135(33) 


4.0427(55) 


-0.8085(11) 


0.8 


7.4920(32) 


13 


2774(32) 


4.8211(54) 


-0.9642(11) 


0.9 


6.6181(33) 


13 


3428(33) 


5.6039(55) 


-1.1208(11) 


1.0 


5.7264(33) 


13 


4109(33) 


6.4038(55) 


-1.2808(11) 


1.1 


4.8031(33) 


13 


4827(33) 


7.2330(55) 


-1.4466(11) 


1.2 


3.8324(34) 


13 


5595(34) 


8.1059(56) 


-1.6212(11) 


1.3 


2.7953(33) 


13 


6428(33) 


9.0396(55) 


-1.8079(11) 


1.4 


1.6679(34) 


13 


7343(34) 


10.0553(56) 


-2.0111(11) 


1.5 


0.4177(33) 


13 


8368(33) 


11.1826(55) 


-2.2365(11) 


1.6 


-1.0018(33) 


13 


9535(33) 


12.4627(55) 


-2.4925(11) 


1.7 


-2.6633(33) 


14 


0898(33) 


13.9609(55) 


-2.7922(11) 


1.8 


-4.6917(33) 


14 


2543(33) 


15.7883(55) 


-3.1577(11) 


1.9 


-7.3491(33) 


14 


4641(33) 


18.1777(55) 


-3.6355(11) 



TABLE IV: Numerical values for Zij (i, j = 1, . . . , 5) as a function of M with plaquette gauge action. 



M 




4^ 








,± 

-'22 


-^23 


^23 


0.1 


-52 


33012(24) 


-39 


2974(68) 


-48 


1580(13) 


-3.0327(66) 


1.0327(66) 


0.2 


-51 


4142(13) 


-39 


9464(53) 


-47 


50286(22) 


-1.4678(66) 


-0.5322(66) 


0.3 


-50 


6609(11) 


-40 


4585(55) 


-46 


960537(34) 


-0.2024(66) 


-1.7976(66) 


0.4 


-50 


0050(14) 


-40 


9157(51) 


-46 


49015(35) 


0.9107(65) 


-2.9107(65) 


0.5 


-49 


41704(64) 


-41 


3522(72) 


-46 


0729(17) 


1.9352(66) 


-3.9352(66) 


0.6 


-48 


8806(11) 


-41 


7879(54) 


-45 


698460(58) 


2.9073(65) 


-4.9073(65) 


0.7 


-48 


3854(18) 


-42 


2366(48) 


-45 


36057(67) 


3.8512(66) 


-5.8512(66) 


0.8 


-47 


92464(38) 


-42 


7100(61) 


-45 


05553(70) 


4.7853(65) 


-6.7853(65) 


0.9 


-47 


49355(14) 


-43 


2182(68) 


-44 


7810(12) 


5.7247(66) 


-7.7247(66) 


1.0 


-47 


0888(16) 


-43 


7734(50) 


-44 


53624(49) 


6.6845(65) 


-8.6845(65) 


1.1 


-46 


707975(79) 


-44 


3876(67) 


-44 


3212(12) 


7.6796(66) 


-9.6796(66) 


1.2 


-46 


3494(12) 


-45 


0765(55) 


-44 


137215(96) 


8.7271(67) 


-10.7271(67) 


1.3 


-46 


011799(83) 


-45 


8593(65) 


-43 


9864(10) 


9.8475(66) 


-11.8475(66) 


1.4 


-45 


6944(11) 


-46 


7608(56) 


-43 


872137(23) 


11.0664(67) 


-13.0664(67) 


1.5 


-45 


3964(14) 


-47 


8155(51) 


-43 


79961(34) 


12.4191(65) 


-14.4191(65) 


1.6 


-45 


11658(35) 


-49 


0719(70) 


-43 


7758(15) 


13.9553(66) 


-15.9553(66) 


1.7 


-44 


8518(17) 


-50 


6049(49) 


-43 


81064(62) 


15.7531(66) 


-17.7531(66) 


1.8 


-44 


59306(72) 


-52 


5390(73) 


-43 


9174(18) 


17.9460(66) 


-19.9460(66) 


1.9 


-44 


3086(13) 


-55 


1218(53) 


-44 


11077(18) 


20.8132(66) 


-22.8132(66) 
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TABLE V: Numerical values for j = 1, . . . , 5) as a function of M with plaquette gauge action. 



M 




■^33 








^44 


7 + 

-'45 


^45 


0.1 


-37 


609(11) 


-40 


2919(79) 


-36 


259(14) 


-0.77960(37) 


1.8980(18) 


0.2 


-39 


3012(97) 


-41 


2018(64) 


-38 


734(13) 


-0.69265(37) 


1.4633(18) 


0.3 


-40 


6569(99) 


-41 


9248(66) 


-40 


722(13) 


-0.62236(37) 


1.1118(18) 


0.4 


-41 


8562(94) 


-42 


5675(61) 


-42 


478(13) 


-0.56052(36) 


0.8026(18) 


0.5 


-42 


976(12) 


-43 


1748(83) 


-44 


110(15) 


-0.50360(36) 


0.5180(18) 


0.6 


-44 


0595(97) 


-43 


7725(65) 


-45 


680(13) 


-0.44959(36) 


0.2480(18) 


0.7 


-45 


1374(92) 


-44 


3785(59) 


-47 


230(13) 


-0.39715(37) 


-0.0142(18) 


0.8 


-46 


234(10) 


-45 


0075(72) 


-48 


793(14) 


-0.34526(36) 


-0.2737(18) 


0.9 


-47 


368(11) 


-45 


6723(79) 


-50 


397(14) 


-0.29307(37) 


-0.5346(18) 


1.0 


-48 


5631(93) 


-46 


3875(61) 


-52 


072(13) 


-0.23975(36) 


-0.8013(18) 


1.1 


-49 


841(11) 


-47 


1675(78) 


-53 


847(14) 


-0.18447(37) 


-1.0777(18) 


1.2 


-51 


228(10) 


-48 


0310(66) 


-55 


758(13) 


-0.12627(37) 


-1.3686(19) 


1.3 


-52 


758(11) 


-49 


0006(76) 


-57 


848(14) 


-0.06403(37) 


-1.6799(18) 


1.4 


-54 


472(10) 


-50 


1052(67) 


-60 


172(13) 


0.00369(37) 


-2.0184(19) 


1.5 


-56 


4282(95) 


-51 


3854(62) 


-62 


804(13) 


0.07884(36) 


-2.3942(18) 


1.6 


-58 


709(11) 


-52 


8977(81) 


-65 


853(15) 


0.16418(37) 


-2.8209(18) 


1.7 


-61 


4402(93) 


-54 


7304(60) 


-69 


483(13) 


0.26406(37) 


-3.3203(18) 


1.8 


-64 


836(12) 


-57 


0300(84) 


-73 


976(15) 


0.38589(37) 


-3.9294(18) 


1.9 


-69 


3305(97) 


-60 


0906(64) 


-79 


904(13) 


0.54518(37) 


-4.7259(18) 



TABLE VL Numerical values for Zij (i, j = 1, . . . , 5) as a function of M with plaquette gauge action. 



M 




^55 


^55 


2 + 

^54 


^54 


0.1 


-50 


2390(38) 


-44.2063(28) 


2.6939(55) 


0.6612(11) 


0.2 


-48 


3667(54) 


-43.8989(12) 


1.3898(55) 


0.9220(11) 


0.3 


-46 


8402(52) 


-43.6378(14) 


0.3353(55) 


1.1329(11) 


0.4 


-45 


5041(54) 


-43.4147(11) 


-0.5922(54) 


1.3184(11) 


0.5 


-44 


2900(34) 


-43.2252(32) 


-1.4460(55) 


1.4892(11) 


0.6 


-43 


1594(51) 


-43.0668(14) 


-2.2561(54) 


1.6512(11) 


0.7 


-42 


0874(58) 


-42.93862(80) 


-3.0427(55) 


1.8085(11) 


0.8 


-41 


0558(43) 


-42.8412(21) 


-3.8211(54) 


1.9642(11) 


0.9 


-40 


0507(39) 


-42.7754(27) 


-4.6039(55) 


2.1208(11) 


1.0 


-39 


0594(56) 


-42.74392(96) 


-5.4038(55) 


2.2808(11) 


1.1 


-38 


0704(40) 


-42.7500(26) 


-6.2330(55) 


2.4466(11) 


1.2 


-37 


0717(53) 


-42.7988(14) 


-7.1059(56) 


2.6212(11) 


1.3 


-36 


0494(41) 


-42.8970(25) 


-8.0396(55) 


2.8079(11) 


1.4 


-34 


9872(52) 


-43.0536(15) 


-9.0553(56) 


3.0111(11) 


1.5 


-33 


8626(54) 


-43.2816(11) 


-10.1826(55) 3.2365(11) 


1.6 


-32 


6439(37) 


-43.5992(29) 


-11.4627(55) 3.4925(11) 


1.7 


-31 


2805(57) 


-44.03354(85) 


-12.9609(55) 3.7922(11) 


1.8 


-29 


6816(33) 


-44.6276(33) 


-14.7883(55) 4.1577(11) 


1.9 


-27 


6450(53) 


-45.4581(13) 


-17.1777(55) 4.6355(11) 
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TABLE VII: Numerical values foruy — 1, ■ • ■ , 5) as a function of M with Iwasaki gauge action. 



M 






^11 








^23 


0.1 


14.1911(73) 


10 


1069(73) 


13 


5104(18) 


4.0842(66) 


-4.0842(66) 


0.2 


13.7193(73) 


11 


0602(73) 


13 


2761(18) 


2.6591(66) 


-2.6591(66) 


0.3 


13.3494(73) 


11 


8105(73) 


13 


0929(18) 


1.5389(66) 


-1.5389(66) 


0.4 


13.0327(72) 


12 


4558(72) 


12 


9365(18) 


0.5769(65) 


-0.5769(65) 


0.5 


12.7480(73) 


13 


0378(73) 


12 


7963(18) 


-0.2899(66) 


0.2899(66) 


0.6 


12.4836(72) 


13 


5804(72) 


12 


6664(18) 


-1.0968(65) 


1.0968(65) 


0.7 


12.2320(73) 


14 


0990(73) 


12 


5431(18) 


-1.8670(66) 


1.8670(66) 


0.8 


11.9874(72) 


14 


6053(72) 


12 


4237(18) 


-2.6179(65) 


2.6179(65) 


0.9 


11.7451(74) 


15 


1082(74) 


12 


3057(18) 


-3.3631(66) 


3.3631(66) 


1.0 


11.5011(73) 


15 


6173(73) 


12 


1871(18) 


-4.1162(65) 


4.1162(65) 


1.1 


11.2508(73) 


16 


1409(73) 


12 


0658(18) 


-4.8901(66) 


4.8901(66) 


1.2 


10.9896(75) 


16 


6892(75) 


11 


9395(19) 


-5.6997(67) 


5.6997(67) 


1.3 


10.7118(73) 


17 


2742(73) 


11 


8055(18) 


-6.5624(66) 


6.5624(66) 


1.4 


10.4101(75) 


17 


9105(75) 


11 


6602(19) 


-7.5004(67) 


7.5004(67) 


1.5 


10.0748(73) 


18 


6191(73) 


11 


4988(18) 


-8.5444(65) 


8.5444(65) 


1.6 


9.6912(74) 


19 


4294(74) 


11 


3142(18) 


-9.7382(66) 


9.7382(66) 


1.7 


9.2365(73) 


20 


3887(73) 


11 


0952(18) 


-11.1523(66) 11.1523(66) 


1.8 


8.6700(73) 


21 


5799(73) 


10 


8217(18) 


-12.9099(66) 12.9099(66) 


1.9 


7.9044(73) 


23 


1801(73) 


10 


4503(18) 


-15.2758(66) 15.2758(66) 



TABLE VIII: Numerical values for mj 1, . . . , 5) as a function of M with Iwasaki gauge action. 



M 








«44 




^45 


0.1 


7.384(12) 


9.4262(84) 


5.342(15) 


0.22690(37) 


-1.1345(18) 


0.2 


9.287(12) 


10 


6170(84) 


7.958(15) 


0.14773(37) 


-0.7386(18) 


0.3 


10.785(12) 


11 


5540(84) 


10.015(15) 


0.08549(37) 


-0.4275(18) 


0.4 


12.071(12) 


12 


3596(83) 


11.783(15) 


0.03205(36) 


-0.1603(18) 


0.5 


13.231(12) 


13 


0862(84) 


13.376(15) 


-0.01610(36) 


0.0805(18) 


0.6 


14.312(12) 


13 


7632(83) 


14.860(15) 


-0.06093(36) 


0.3047(18) 


0.7 


15.344(12) 


14 


4102(84) 


16.277(15) 


-0.10372(37) 


0.5186(18) 


0.8 


16.351(11) 


15 


0416(83) 


17.660(15) 


-0.14544(36) 


0.7272(18) 


0.9 


17.350(12) 


15 


6687(85) 


19.032(15) 


-0.18684(37) 


0.9342(18) 


1.0 


18.361(12) 


16 


3033(84) 


20.420(15) 


-0.22868(36) 


1.1434(18) 


1.1 


19.401(12) 


16 


9559(84) 


21.846(15) 


-0.27167(37) 


1.3584(18) 


1.2 


20.489(12) 


17 


6392(86) 


23.339(15) 


-0.31665(37) 


1.5832(19) 


1.3 


21.649(12) 


18 


3680(84) 


24.930(15) 


-0.36458(37) 


1.8229(18) 


1.4 


22.911(12) 


19 


1606(86) 


26.661(15) 


-0.41669(37) 


2.0834(19) 


1.5 


24.315(12) 


20 


0432(84) 


28.588(15) 


-0.47469(36) 


2.3734(18) 


1.6 


25.921(12) 


21 


0524(85) 


30.791(15) 


-0.54101(37) 


2.7051(18) 


1.7 


27.824(12) 


22 


2475(84) 


33.400(15) 


-0.61957(37) 


3.0979(18) 


1.8 


30.186(12) 


23 


7315(84) 


36.641(15) 


-0.71721(37) 


3.5861(18) 


1.9 


33.364(12) 


25 


7261(84) 


41.002(15) 


-0.84865(37) 


4.2433(18) 
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TABLE IX: Numerical values for Vij (i, j = 1, . . . , 5) as a function of M with Iwasaki gauge action. 



M 


4, 




"55 


^^54 


"54 


0.1 


16.6870(33) 


12 


6028(33) 


-3.4035(55) 


0.6807(11) 


0.2 


15.3443(33) 


12 


6852(33) 


-2.2159(55) 


0.4432(11) 


0.3 


14.2899(33) 


12 


7510(33) 


-1.2824(55) 


0.2565(11) 


0.4 


13.3852(32) 


12 


8083(33) 


-0.4808(54) 


0.0962(11) 


0.5 


12.5708(33) 


12 


8607(33) 


0.2416(55) 


-0.0483(11) 


0.6 


11.8133(33) 


12 


9102(33) 


0.9140(54) 


-0.1828(11) 


0.7 


11.0910(33) 


12 


9580(33) 


1.5559(55) 


-0.3112(11) 


0.8 


10.3875(32) 


13 


0055(32) 


2.1816(54) 


-0.4363(11) 


0.9 


9.6899(33) 


13 


0530(33) 


2.8025(55) 


-0.5605(11) 


1.0 


8.9856(33) 


13 


1018(33) 


3.4302(55) 


-0.6860(11) 


1.1 


8.2624(33) 


13 


1525(33) 


4.0751(55) 


-0.8150(11) 


1.2 


7.5064(34) 


13 


2061(34) 


4.7497(56) 


-0.9499(11) 


1.3 


6.7014(33) 


13 


2638(33) 


5.4687(55) 


-1.0937(11) 


1.4 


5.8266(34) 


13 


3270(34) 


6.2503(56) 


-1.2501(11) 


1.5 


4.8532(33) 


13 


3976(33) 


7.1203(55) 


-1.4241(11) 


1.6 


3.7400(33) 


13 


4782(33) 


8.1152(55) 


-1.6230(11) 


1.7 


2.4212(33) 


13 


5735(33) 


9.2936(55) 


-1.8587(11) 


1.8 


0.7807(33) 


13 


6905(33) 


10.7582(55) 


-2.1516(11) 


1.9 


-1.4308(33) 


13 


8449(33) 


12.7298(55) 


-2.5460(11) 



TABLE X: Numerical values for Zij {i, j = 1, ... ,5) as a function of M with Iwasaki gauge action. 



M 




4i 




Zll 




^22 


2 + 
^23 


^23 


0.1 


-30 


25980(58) 


-15 


1756(72) 


-25 


7458(17) 


-5.0842(66) 


3.0842(66) 


0.2 


-29 


43439(58) 


-15 


7753(72) 


-25 


1579(17) 


-3.6591(66) 


1.6591(66) 


0.3 


-28 


77133(70) 


-16 


2324(73) 


-24 


6815(18) 


-2.5389(66) 


0.5389(66) 


0.4 


-28 


20508(57) 


-16 


6282(71) 


-24 


2756(17) 


-1.5769(65) 


-0.4231(65) 


0.5 


-27 


70614(68) 


-16 


9960(72) 


-23 


9211(18) 


-0.7101(66) 


-1.2899(66) 


0.6 


-27 


25798(48) 


-17 


3548(70) 


-23 


6075(16) 


0.0968(65) 


-2.0968(65) 


0.7 


-26 


85022(60) 


-17 


7173(72) 


-23 


3281(17) 


0.8670(66) 


-2.8670(66) 


0.8 


-26 


47585(51) 


-18 


0938(70) 


-23 


0788(16) 


1.6179(65) 


-3.6179(65) 


0.9 


-26 


12985(57) 


-18 


4929(72) 


-22 


8570(17) 


2.3631(66) 


-4.3631(66) 


1.0 


-25 


80864(62) 


-18 


9249(72) 


-22 


6613(17) 


3.1162(65) 


-5.1162(65) 


1.1 


-25 


50947(72) 


-19 


3996(73) 


-22 


4912(18) 


3.8901(66) 


-5.8901(66) 


1.2 


-25 


23026(60) 


-19 


9299(73) 


-22 


3469(17) 


4.6997(67) 


-6.6997(67) 


1.3 


-24 


96938(47) 


-20 


5318(71) 


-22 


2298(16) 


5.5624(66) 


-7.5624(66) 


1.4 


-24 


72536(72) 


-21 


2258(74) 


-22 


1421(18) 


6.5004(67) 


-8.5004(67) 


1.5 


-24 


49673(63) 


-22 


0411(72) 


-22 


0875(17) 


7.5444(65) 


-9.5444(65) 


1.6 


-24 


28127(49) 


-23 


0195(71) 


-22 


0710(16) 


8.7382(66) 


-10.7382(66) 


1.7 


-24 


07481(51) 


-24 


2271(71) 


-22 


1002(16) 


10.1523(66) 


-12.1523(66) 


1.8 


-23 


86686(71) 


-25 


7767(73) 


-22 


1852(18) 


11.9099(66) 


-13.9099(66) 


1.9 


-23 


62377(52) 


-27 


8995(71) 


-22 


3364(16) 


14.2758(66) 


-16.2758(66) 
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TABLE XI; Numerical values for Zij {i, j = 1, . . . , 5) as a function of M with Iwasaki gauge action. 



M 




^33 




7 + 


^44 


z+ 
^45 


^45 


0.1 


-12 


119(12) 


-15 


8282(82) 


-9.744(15) 


-0.89357(37) 


2.4678(18) 


0.2 


-13 


669(12) 


-16 


6654(83) 


-12.006(15) 


-0.81440(37) 


2.0720(18) 


0.3 


-14 


873(12) 


-17 


3093(84) 


-13.770(15) 


-0.75216(37) 


1.7608(18) 


0.4 


-15 


910(11) 


-17 


8654(82) 


-15.288(15) 


-0.69872(36) 


1.4936(18) 


0.5 


-16 


856(12) 


-18 


3777(83) 


-16.668(15) 


-0.65056(36) 


1.2528(18) 


0.6 


-17 


753(11) 


-18 


8709(81) 


-17.968(15) 


-0.60573(36) 


1.0287(18) 


0.7 


-18 


629(12) 


-19 


3618(83) 


-19.229(15) 


-0.56294(37) 


0.8147(18) 


0.8 


-19 


506(11) 


-19 


8634(81) 


-20.481(15) 


-0.52123(36) 


0.6061(18) 


0.9 


-20 


402(12) 


-20 


3867(83) 


-21.750(15) 


-0.47983(37) 


0.3992(18) 


1.0 


-21 


336(12) 


-20 


9442(83) 


-23.060(15) 


-0.43799(36) 


0.1899(18) 


1.1 


-22 


326(12) 


-21 


5479(84) 


-24.438(15) 


-0.39499(37) 


-0.0250(18) 


1.2 


-23 


396(12) 


-22 


2132(84) 


-25.913(15) 


-0.35002(37) 


-0.2499(19) 


1.3 


-24 


573(11) 


-22 


9589(82) 


-27.521(15) 


-0.30209(37) 


-0.4896(18) 


1.4 


-25 


893(12) 


-23 


8092(85) 


-29.310(15) 


-0.24998(37) 


-0.7501(19) 


1.5 


-27 


404(12) 


-24 


7985(83) 


-31.343(15) 


-0.19198(36) 


-1.0401(18) 


1.6 


-29 


178(12) 


-25 


9758(82) 


-33.714(15) 


-0.12566(37) 


-1.3717(18) 


1.7 


-31 


329(11) 


-27 


4191(82) 


-36.571(15) 


-0.04710(37) 


-1.7645(18) 


1.8 


-34 


050(12) 


-29 


2617(84) 


-40.172(15) 


0.05055(37) 


-2.2527(18) 


1.9 


-37 


750(11) 


-31 


7788(82) 


-45.055(15) 


0.18199(37) 


-2.9099(18) 



TABLE XIL Numerical values for z^j (i, j = 1, . . . , 5) as a function of M with Iwasaki gauge action. 



M 




4 




255 


z+ 
-^54 


^54 


0.1 


-29 


4224(34) 


-21 


3382(31) 


4.4035(55) 


0.3193(11) 


0.2 


-27 


7261(34) 


-21 


0670(31) 


3.2159(55) 


0.5568(11) 


0.3 


-26 


3784(33) 


-20 


8395(33) 


2.2824(55) 


0.7435(11) 


0.4 


-25 


2243(34) 


-20 


6474(31) 


1.4808(54) 


0.9038(11) 


0.5 


-24 


1957(33) 


-20 


4855(32) 


0.7584(55) 


1.0483(11) 


0.6 


-23 


2544(35) 


-20 


3512(30) 


0.0860(54) 


1.1828(11) 


0.7 


-22 


3759(34) 


-20 


2430(32) 


-0.5559(55) 


1.3112(11) 


0.8 


-21 


5427(34) 


-20 


1606(30) 


-1.1816(54) 


1.4363(11) 


0.9 


-20 


7413(35) 


-20 


1044(31) 


-1.8025(55) 


1.5605(11) 


1.0 


-19 


9598(34) 


-20 


0761(32) 


-2.4302(55) 


1.6860(11) 


1.1 


-19 


1877(33) 


-20 


0778(33) 


-3.0751(55) 


1.8150(11) 


1.2 


-18 


4138(35) 


-20 


1135(32) 


-3.7497(56) 


1.9499(11) 


1.3 


-17 


6257(36) 


-20 


1881(30) 


-4.4687(55) 


2.0937(11) 


1.4 


-16 


8085(34) 


-20 


3088(33) 


-5.2503(56) 


2.2501(11) 


1.5 


-15 


9418(34) 


-20 


4862(32) 


-6.1203(55) 


2.4241(11) 


1.6 


-14 


9968(36) 


-20 


7350(31) 


-7.1152(55) 


2.6230(11) 


1.7 


-13 


9262(35) 


-21 


0785(31) 


-8.2936(55) 


2.8587(11) 


1.8 


-12 


6442(33) 


-21 


5540(33) 


-9.7582(55) 


3.1516(11) 


1.9 


-10 


9552(35) 


-22 


2310(31) 


-11.7298(55) 3.5460(11) 



20 



TABLE XIII: Numerical values of (i, j = 1, . . . , 5) for plaquette gauge action (upper part) and Iwasaki 
gauge action (lower part). 



6.0 



M 

1.3112 



P 

0.59374 



2.1793 1.80 



7+ 
^11 

0.7287 

7+ 

0.6325 

^4"5 



0.7289 

^32 

±0.0207 



-0.02371 0.8674 



^22 

0.7564 
0.6853 
0.7710 



7+ 
^23 

0.1378 

^44 

0.5613 

7+ 

-0.1125 



^23 

-0.1654 

^4+5 

-0.0007790 

^54 

0.03906 



/3 

2.6 



M 

1.4198 



P 

0.67063 



R 

0.45283 



2.2479 1.80 



0.8062 

7+ 
-^33 



0.8531 



^32 



^22 

0.8425 
7+ 



^23 

0.09548 

^44 



0.7847 ±0.02135 0.8158 0.7346 



^45 



^55 



7+ 



^23 

-0.1239 
7+ 

-0.003395 



^54 



-0.01150 0.9207 0.8680 -0.07719 0.03252 



21 



a,i 




